function msvalue(num)
 load('orthbase');
% 
% load('betas');
% 
% row=size(betas,1); colum=size(betas,2);
% if row >=2
%     rest(row)=0;
% end
% tmp = JNKR(sum(betas,2)==1,:);

time=linspace(0,0.5,num+1);

msval=zeros(num+1,100);
for ii=1:num+1
    for jj=1:100
        msval(ii,jj)=ms{jj}(time(ii));
    end
end

fname='msval';
save(fname,'msval');
